On Quantum Correlations across the Black Hole Horizon 
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Inspired by the condensed-matter analogues of black holes, we study the quantum correlations 
across the event horizon reflecting the entanglement between the outgoing particles of the Hawking 
radiation and their in-falling partners. For a perfectly covariant theory, the total correlation is 
conserved in time and piles up arbitrary close to the horizon in the past, where it merges into 
the singularity of the vacuum two-point function at the light cone. After modifying the dispersion 
relation (i.e., breaking Lorentz invariance) for large k, on the other hand, the light cone is smeared 
out and the entanglement is not conserved but actually created in a given rate per unit time. 

PACS numbers: 04.62. +v, 04.70.Dy, 04.60.-m. 



I. INTRODUCTION 

One of the remaining mysteries of modern physics is 
the question of why black holes seem to behave as ther- 
mal objects, described by the Hawking temperature 
The thermal nature of black holes entails important con- 
cepts such as black hole entropy [2] and the information 
"paradox" . It is now widely believed that understanding 
the origin of thcrmality would be a major step towards 
unifying gravity and quantum theory. 

Particularly puzzling is the fact that this thermal na- 
ture is apparently not caused by some sort of equilibra- 
tion process, but by a dynamical quantum mechanism. 
In the semi-classical description, the state of the quan- 
tum fields propagating in the gravitational background 
is still a pure state. The thermal nature of the Hawk- 
ing radiation is explained by the quantum correlations 
(across the horizon) between the Hawking particles and 
their infalling partners 

For Hawking radiation itself, the origin of the created 
particles and its robustness against modifications of the 
microscopic structure have been studied in many pub- 
lications - often inspired by the condensed-matter ana- 
logues of black holes ("dumb holes"), see, e.g., UHl]- In 
contrast, the correlations across the horizon, their origin 
and dependence on the microscopic structure have been 
studied in far less detail mi- 



ll. TWO-POINT FUNCTION 

Exploiting the fact that Hawking radiation is basically 
a one-dimensional effect and applies to all fields, we study 
a massless scalar field cf> in a 1+1 dimensional space-time 
described by the Painleve-Gullstrand-Lemaitre coordi- 
nates (h = C = Gn = fcj3 = 1) 



ds 2 



dt 2 
[1- 



- [dx + X>(x)dt} 2 

V 2 (x)]dt 2 - 2v(x)dtdx - dx 2 



where t)(x) can be visualized as the local velocity of 
freely falling frames as measured from infinity. The 
Schwarzschild metric is obtained by = \j2Mjr, but 
we shall consider arbitrary profiles t>(x). In the standard 
manner, we introduce light-cone variables u and v via 



ds 2 = [1 - t> 2 (x)] [dt - 
= [1 - 2 (x)] dudv . 



dx 



1 - t>(x) 



dt 



dx 



1 + o(x) 



(2) 



with u diverging at the (future) event horizon = 1. 
The past horizon with = — 1 corresponds to v f oo. 
For future convenience, we introduce the tortoise coor- 
dinate dx* = dx/[l — t>(x)] with u = t — x*. After the 
standard transformation to regular coordinates for < 1 
(i.e., outside the horizon) 



U 



1 



V = 



1 



(3) 



where k = (dv/dx) horizon i s the surface gravity, and ana- 
lytic continuation beyond horizon vj > 1 where U > 0, we 
obtain the line element with a regular conformal factor 



ds 2 =l3 2 cgu UUV)dUdV. 



(4) 



Due to the conformal invariance of the massless scalar 
field in 1+1 dimensions, we may directly read off the 
two-point function(s). In the Boulware state [8], which 
is the ground state of the Hamiltonian generating the t- 
evolution (i.e., of all stationary observers), it behaves as 
ln(AwAu). However, this quantity is clearly divergent at 
both horizons, u j oo and v t oo. Black hole evaporation 
is described by the Unruh state, which is regular across 
the black hole (future) horizon at U — leading to the 
two-point function oc In(AUAv). In this case, the ingoing 
u-modes are in their ground state. The Israel-Hartle- 
Hawking state [t| is regular across both horizons U = 
and V = and is thermal in both directions. Its two- 
point function reads (up to an undetermined constant 
reflecting the IR-divergence in 1+1 dimensions) 



(1) 



(<j>{JJ,V)<l>{U' ,V')) 



1 

47T 



ln(AZ7AV) . 



(5) 
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where AU = U — U' and AV — V — V. Since the ingoing 
(V, v) sector is decoupled from the outgoing (U, u) sector 
(in this 1+1 dimensional set-up), we focus on the relevant 
AU part (•••)[/ describing the Hawking radiation and the 
in-falling partner particles in the following. If U and U' 
lie on different sides of the horizon UU' < 0, we obtain 



(<f>(U,V)<l>(U' ) V')) u = - 



47T V 



(6) 



In this form, correlations across the horizon do not be- 
come particularly apparent, but if we calculate, for ex- 
ample, the correlator 



^(t,xW,x'))u = -^d t d t , m(e- 
k 2 1 



16tt cosh 2 (K[u + u']/2) ' 



(7) 



with u + u' = t + 1' — — x' r , we see that it has a peak if 
we regard the the Hawking particle at (x*) at time t and 
its in- falling partner is at (x[) at time t' Q- 

Due to u f oo at the horizon, this correlator vanishes 
where = 1. This reflects the critical slow-down (of 
the u-modes) at the horizon in terms of the ^coordinate. 
However, for other quantities such as the momentum den- 
sity II = {pt — od x )cf), the correlator does not vanish when 
approaching the horizon 



(n(t 1 x)n(i',a; , ))i/ = -^d x d x ,ln(e- KU + e- KU ') 
K 2 [l-0(a;)]- 1 [l-0(x')]- 1 



16tt cosh 2 (k[w + u']/2) 



(8) 



in view of c^it = — 1/(1 — t)). We observe that the (ft cor- 
relator across the horizon is always positive whereas the 
n correlator is negative. This can be explained by the 
fact that the particles of the Hawking radiation and their 
in-falling partners have the same conserved frequency to 
measured with respect to the time t but their fc-values 
have the opposite sign. Note that ui and k should not 
be confused with energy and momentum: The Hawking 
particle has positive energy and momentum whereas the 
in-falling partner has negative energy but positive mo- 
mentum, see, e.g., [Toj . 



III. INTEGRALS 

Since the n-corrclator is total derivative, we may easily 
derive the total correlation integrated from the horizon 
at x — up to spatial infinity 



dx(U(t,x)U(t',x' < 0)) t 



l 



4tt 1 - o(x') ' 



(9) 



and we find that it is independent of t — t'. Consequently, 
for each point x' < inside the horizon (at x = 0), the in- 
tegral of the correlations across the horizon is conserved, 



i.e., independent of the time-slice. In the far future 1 1 oo, 
the Hawking particles carry the correlations to spatial 
infinity - as one would expect. In the past t I — oo, 
however, the correlations pile up near the horizon where 
u ~ t — k^ 1 Iti(kx) 



(IL(t I -oo, x I 0)n(t', x' < 0))u « 

K [l-P^')]- 1 1 

167T cosh 2 (At[u + u']/2) x 



(10) 



i.e., they are concentrated in a small spatial volume and 
have a large amplitude (in order to keep the integral con- 
stant). As we shall see in the next Section, these piled- up 
correlations merge into the singularity of the two-point 
function at the light cone (which approaches the horizon) 
and become virtually indistinguishable from the quantum 
vacuum fluctuations. For a modified dispersion relation, 
this picture changes drastically, see Section IVl 
Similarly, we may evaluate the time-integral 



In contrast to the x-integral above, this result will survive 
for a modified dispersion relation (apart from some cor- 
rections at short length scales, see Section|V|. Of course, 
we may also derive the time-integral for 



+ oo 



dt(<P(t,xW,x / )) u = — 

47T 



(12) 



The difference between the two expressions |TT|) and (fT"2")) 
can again be traced back to the fact that the Hawking 
particles and their in- falling partners have the same to, 
but different k depending on their positions x and x' 
(gravitational red-shift u = [1 — t)]fc). 



IV. ANALYTICALLY SOLVABLE EXAMPLE 

It might be illustrative to apply the above formulae 
(which are valid for arbitrary profiles u) to some simple 
example which allows us to write down closed expres- 
sions. To this end, let us choose the velocity profile 



0(x) 



1 



7 



• tanh(7ir) . 



(13) 



where 7 is some parameter and k the surface gravity. In 
this case, the light cone coordinates read 

u = t - - ln[sinh(7ar)l U = --e~ Kt sinh(7a;) , (14) 

K K 

and the correlator for all x and x' becomes 



{4>(t,x)4>(t,xr)) u = -— u 



k 2 sinh(7a;) sinh(7o:') 



4-7T [sinh(7a;) — sinh(7x')] 



(15) 
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This function is positive on opposite sides of the horizon 
xx' < 0, negative when x and x' lie on the same side 
xx 1 > 0, and goes to zero if one of the two points ap- 
proaches the horizon at x — (critical slow-down). It 
reproduces the usual l/(x — x 1 ) 2 singularity at x = x' 
and vanishes asymptotically x — > ±oo and x' — > ±oo. 
The correlations across the horizon manifest themselves 
in the global maximum at x = —x' with t = t' . 

Although one gets the same l/(x — x') 2 singularity at 
x = x' and the same asymptotic behavior for x — > ±oo 
and x' — > ±oo, the situation is a bit different for the 
canonical momentum density II = — 00' 



(U{t,x)U{t,x'))u = 



7 2 cosh^x) cosh(7x') 
47r [sinh(7x) — sinh^a;')] 2 



.(16) 



First, this correlator is negative everywhere and does not 
vanish at the horizon. Second, we get a local minimum 
at x = —x' only if \x\ > x* rs 1-8/7, i- e -, f &r enough 
away from the horizon. Close to the horizon, the local 
minimum merges into the light-cone singularity and dis- 
appears, see Figure [TJ This absence of structure near the 
horizon can be interpreted as further confirmation of the 
picture that the Hawking radiation is not created very 
close to the horizon but rather in a region of finite spa- 
tial extent 0(1/7). Note that this length scale I/7 could 
in principle be quite different from the scale set by the 
Hawking temperature, i.e., the surface gravity k. 

It should also be mentioned that the (^-correlator does 
have a maximum arbitrarily close to the horizon, but the 
maximum value becomes very small. Consistent with this 
observation, the correlation conservation law (j9|) applies 
to the n-correlator - for the ^-correlator, one would need 
an additional integrating factor 1/[1 — o(cc)]. 



V. MODIFIED DISPERSION 

In the previous Sections, we found the infinite pile up 
of correlation close to the horizon in the past, such that 
the correlations constantly emerge out of the singularity 
at the light-cone (which approaches the null surface of the 
horizon in the past). However, modifying the dispersion 
relation at large k (e.g., motivated by the condensed- 
matter analogues of black holes 041]); we expect this 
behavior to change: First, a modified dispersion relation 
smears out the light cone such that the two-point func- 
tion typically becomes regular everywhere except at the 
space-time coincidence point x = x' and t = t'. Second, 
tracing the particles of the Hawking radiation back in 
time, they do not originate from an small vicinity of the 
horizon in the presence of a modified dispersion relation 

M- 

In order to deal with a solvable example, let us switch 
to the Eddington-Finkelstein coordinates (v,r) 



ds 2 



dv — 2dv dr . 



r I 



(17) 




FIG. 1: Plot of the correlator (| ltjp as a function of position x 
outside the horizon for various values of the partner position 
x' inside the horizon (both in units of 7). For large values of 
x (i.e., far from the horizon) the correlation as a function of x 
is simply a bump of width about 2 and constant height which 
is located at x — x' . For small x' it merges into the general 
vacuum fluctuations near the horizon, i.e., the l/(x — x') 2 
singularity. I.e., the correlation of the "Hawking" particles 
with their infalling partners becomes visible only if the two 
points are well separated from the horizon, even though the 
integral (J5J) of the correlation over a; is a constant. 



and modify the dispersion relation by inserting a general 
function f(d 2 ) containing higher-order spatial derivatives 
into the corresponding action for a scalar field 



C = -(d v 4>)d r 



1 



2M 



f(d 2 r 



d r 4>. (18) 



This action allows us to define a conserved inner product 
(0i \<f> 2 ) = i [ <t>$ d» 02 = % [ dr cf>l d r 2 , (19) 



where cf>l dr 02 = 0i<9 r 02 — 02<9 r 0*. The momentum 
density n = — d r (j) satisfies the wave equation 



2d v + d r 



r 



n = 0. 



(20) 



which is of first order in d v and hence automatically se- 
lects the outgoing u-sector only. In a stationary state, 
the two-point function can be Fourier expanded via 

{n(v, r)lKV, r')> = / du e-^ v ~ v "> 9ul {r, r') , (21) 
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where gu>(r,r') solves the ordinary differential equation 



-2iu + d r 



2M 
r 



f(d 2 ) 



gUr,r') = Q (22) 



for r and the same for r' with +2iuj. Assuming that this 
differential equation together with the asymptotic condi- 
tions (freely falling ground state for large k at all posi- 
tions r) uniquely determines the r-dependence of g u (r, r') 
(and thus the same for r'), we find that g w (r, r') factorizes 



9u(ry) = h u (r)hZ(r'). 



(23) 



Here, we are interested in the region near the horizon 
(where the pile-up of correlation occurred) and thus we 
employ the near-horizon approximation 



1-™=2kx+ 0(k 2 x 2 ), 
r 



(24) 



resulting in the (approximate) differential equation which 
can be solved via a Fourier-Laplace transformation 



(-2ioj + d x [2nx + f{dl)})K{x) = 
(-2iuj + ik[2nid k + f{-k 2 )])h u (k) = 0. 



(25) 



This is now a first-order ordinary differential equation in 
k and its general solution can be written as 



K{x) = dk k^l K exp \ ikx - -F{k 2 ) 



(26) 



different choices for the branch cut from fc~'"/ K . If the 
branch cut lies in the upper complex half plane 9(fc) > 0, 
the solution connects the final Hawking mode to the ini- 
tial positive/negative (pseudo) norm modes @ 

a^t (x) + Pu€ (x) -> <C wldns (z > 0) , (27) 

where (x) has positive (pseudo) norm (|19|) and 
4>u (x) negative (pseudo) norm. On the other hand, the 
branch cut in the lower complex half plane 3(/c) < 
connects the final mode of the in-falling partners to the 
initial positive/negative (pseudo) norm modes 



partner 



(a; < o) 



(28) 



Assuming that the quantum state we have corresponds to 
the freely falling ground state for large k 7 the asymptotic 
condition for hu,{x") implies that it has no contribution 
from 0k (x). Thus, we take a suitable linear combination 
of the two solutions (j27|) and (|28| with the branch cut 
in the upper and lower complex half-plane, respectively, 
which yields 



h u (x) = A/L / dk k~ iu/n exp \ ikx - -F(k) 



(29) 



where AC is a normalization factor. This expression is 
quite natural since the boundary condition (freely falling 
ground state for large k) implies that h u (x) only con- 
tains positive fc-values with a positive (pseudo) norm 
(fT9l) . Coarse-graining over large length scales, we do 
not see the impact of F(k) and this function behaves as 
h(x) ~ xl 1 ^/" -1 , but on short distances, it also contains 
the rapidly oscillating in- mode <f>* (x). 



where £ is an appropriate contour in the complex plane 
and dF/dk — f(—k 2 )/2 accounts for the modified dis- 
persion relation. In order to determine the correct inte- 
gration contour in the complex plane, we have to study 

I 

Now we are in the position to study the full correlator. Inserting ([29| into (|2Tj) and (|23|) . the total expression reads 

+ oo oo oo 

k' 



(n(v,x)U(v',x')) 



due-^ v - v ' ] LA/?! / dk / dk' exp 



i— In — 

k k 



+ ikx — ik x 



-F(k) + -F(k') 

K K 



(30) 



The cj-integral yields the Fourier transform of \Af 2 \, which we denote by A/". Finally, introducing the new variable 
X = e~ K< - v ~ v ^k/k' ' , the integrated correlation across the horizon in analogy to © yields 



dz<n(« ) a;)ri(t; / ,:fi / < 0)) = / dk' exp I -ik'x' + ^F(k') \ x / ^ Af(ln X ) exp | -^F (xfc'e" re(W) ) \ ■ (31) 



For a given point x' < inside the horizon, Eq. (|3"Tj) 
yields the total correlation between that point x' < and 
all positions x outside the horizon up to spatial infinity. 
Setting F — 0, we rederive the result © for o(x) = kx. 
In the far future (v — v') t oo, the exponential pre- factor 



e -n(v-v ) j n -gq^ (|3l|) vanishes and thus the x-integral be- 
comes independent of k' . The remaining /c'-integral then 
just yields h^ =0 (x'). After coarse-graining over large 
length scales, this scales as 1/x'. As a result, we get ba- 
sically the same conservation law as in ([§]) . This is quite 
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natural since it just reflects that fact that the Hawking 
particles carry the correlation away to spatial infinity. 

In the far past (v — v') 4- — oo, on the other hand, the 
exponential pre-factor e - K -( v ~ v ) diverges and thus the 
//-integral in Eq. (f3l"j) is exponentially suppressed due to 
the rapidly oscillating phase F(xk'e~ K ^ v ~ v '). This can 
most easily be seen by changing the integration variable 
to k" = k'e" K, ^ v ~ v >. In the resulting double integral over 
k" and x> the first exponent in Eq. (|3"Tj) can be neglected 
and the integral scales as e K ^ v ~ v \ Consequently, in con- 
trast to the perfectly covariant case (which implies an 
unbounded red-shift near the horizon) , the total correla- 
tion is not conserved in this case (i.e., it does not pile up 
at the horizon) but created at finite times. 

Note that the integral over the full cc-interval (i.e, 
from — oo to +00) vanishes, since II is a spatial deriva- 
tive. Even with a modified dispersion relation, the time- 
integral factorizes exactly 

+00 

J dv{R(v,x)U(v',x')}=27rh u=0 (x)hZ =0 {x'), (32) 
—00 

and, after coarse-graining over large length scales, it be- 
haves as l/\xx'\. 



away from the horizon. For black hole analogues ( "dumb 
holes") in Bose-Einstein condensates, the scalar field 4> 
reflects the phase fluctuations while the momentum den- 
sity II corresponds to the density fluctuations 5g, cf. Q- 

For the II correlator, we found a conservation law (|9]) 
for the total correlation across horizon in a perfectly co- 
variant theory (up to arbitrarily small length scales). 
This means that the correlation to be carried away by 
the Hawking particles in the future must pile up arbi- 
trarily close to the horizon in the past. However, below a 
minimum length scale set by the geometry (not necessar- 
ily the surface gravity) , this piled-up correlation becomes 
virtually indistinguishable from the vacuum singularity 
of the two-point function at the light cone, see Figure [T] 

After modifying the microscopic structure via intro- 
ducing a non-linear dispersion relation at short distances, 
this picture changes drastically: In this case, the corre- 
lation carried away by Hawking radiation in the future 
cannot be traced back to arbitrarily early times and a 
small vicinity of the horizon. As a result, the entangle- 
ment is not conserved but actually created dynamically 
at a finite time. These findings further support the view 
that Hawking radiation is not created at arbitrarily small 
length scales but at finite distances and could be relevant 
for the black hole information "paradox" etc. 



VI. CONCLUSIONS 

We have studied the evolution of the quantum corre- 
lations across the black hole horizon. Since the quantum 
state under consideration is a pure state, all correlations 
imply entanglement and thus entanglement entropy etc. 
Both, the <j> and the n correlator across the horizon ([7]) 
and p]). possess a peak at u = —v! if we are far enough 
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